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Abstract 

Record statistics is the study of how new highs or lows are created and sus¬ 
tained in any dynamical process. The study of the highest or lowest records 
constitute the study of extreme values. This paper represents an exploration 
of record statistics for certain aspects of the classical and quantum stan¬ 
dard map. For instance the momentum square or energy records is shown 
to behave like that of records in random walks when the classical standard 
map is in a regime of hard chaos. However different power laws is observed 
for the mixed phase space regimes. The presence of accelerator modes are 
well-known to create anomalous diffusion and we notice here that the record 
statistics is very sensitive to their presence. We also discuss records in ran¬ 
dom vectors and use it to analyze the quantum standard map via records in 
their eigenfunction intensities, reviewing some recent results along the way. 


1. Introduction 

Breaking a record, or setting a new one has been a human passion, per¬ 
haps one may say weakness, for a while now. Evidence how one yearns for 
a record to be set by a particular sportsperson or the euphoria (depression) 
that sets in when the stock market fluctuations reach a maximum (mini¬ 
mum) never seen before. Of course our view of nature also abounds in such 
characterizations: the highest mountains, the deepest oceans, the smallest 
temperatures reached and so on. While these may seem to be questions 
about extremes, records are a record of extremes, the times they last and the 
number of new ones created. For engineers designing a dam, one important 
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piece of information is when and by how much the water level of the river 
has exceeded all its previous values; similar questions are naturally asked 
about rainfall when planning for agricultural policies is taken up. Practi¬ 
cal considerations such as these have been the early applications in a study 
of events which exceed themselves in some aspect and are often called as 
“record statistics”, akin to the study of extreme statistics in many ways. An 
accessible introduction to record statistics is found in iSchmittma/nn and Zial 
fll999l ). 

Questions about how the records increase with time, or the number of 
records set, are of natural interest in all these complex, sometimes social 
contexts, have therefore been studied for example, Gembris et al. (2002); 


Vogel et al. (2001); Redner and Petersen (2006). A mathematical theory 
of records for independent identically distributed ( i.i.d .) random variables 
has been developed since the pioneering work of Renyi (1962) for example 
developed in Glick (1978); Arnold et al. (1998). The applications in Physics 
started somewhat later, but have by now found use in various problems 
related to random walks, spin-glasses, type II superconductors and quantum 
chaos. In this work, our emphasis is on the record statistics in deterministic 
dynamical systems, both classical and quantum. In particular we choose the 
most studied low-dimensional paradigm of Hamiltonian chaos, namely the 
standard map, or known variously as the kicked-rotor, Chirikov-Taylor map 
etc.. Its quantization has also been extensively studied, as well as realized 
experimentally in cold atom set ups in Moore et al. (1995). 

Now to define records more precisely, given that {X t ,t = 1, • • • , N} is a 
finite time series, the first element, R(l), of the corresponding records series 
is X\ itself and at subsequent times t it will be R(t) = max(A" t , R(t— 1)). As 
X t is a random variable, so is R(t) and properties of this random variable is 
of interest. Thus R(t) which consists of elements of X t in a non-decreasing 
order are the upper record sequence. If the minimum value is taken instead of 
the maximum the record sequence is called the lower record sequence. From 
the definition of records, it is clear that last upper (lower) record is also the 
global maximum (minimum) of the sequence {Xj}. Hence, the statistics of 
last record will correspond to similar results of maximum or minimum from 
extreme value theory. A word of caution is warranted here, as the second 
last upper (lower) record will not be the second maximum (minimum) of the 
sequence. This can be easily understood by visualizing a sequence X t whose 
global maximum (minimum) occurs before second maximum (minimum). As 
in that case, record statistics will not sense the presence of second maximum 
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(minimum). To develop more familiarity with record sequence, let’s take 
sinx, 0 < x < 2tt as an example. In the range [0, 7t/2] the upper record 
sequence is sin x itself and beyond that the entries will be constant 1. 

Our motivations in this study are many. Firstly, while there is an elegant 
and simple theory of records for the case of independent random variables, 
we can naturally expect interesting departures for correlated processes. One 
such class of problems, namely random walks, have been studied in this 
context and remarkable departures and universalities have been uncovered. 
One of our motivations is to see how much of this survives for deterministic 
dynamical systems. For instance it has been appreciated now for more than 
40 years that deterministic systems such as the kicked rotor can display 
normal diffusion, as well as anomalous diffusion. It is then of interest to 
ask if the record statistics behaves in the same manner as that of random 
walks in this scenario, and what happens in the different diffusion regimes. 
Also related question is how do the record statistics change as the system 
undergoes a transition to chaos. The latter question is natural here, but 
not so in the setting of random walks. We also wish to go on to explore 
quantum dynamics and review some of our recent works on record statistics. 
In this publication, all the material pertaining to record statistics in the 
classical standard map or kicked rotor is new, so also detailed derivations of 
the record statistics for delta correlated variables, as well the lower record 
statistics. Some other material, having to do with eigenfunction statistics 
have been partially presented in publications before (see 
(2013.)) and is included here as a brief review. 


Srivastava et al. 


1.1. Independent and identically distributed (i.i.d.) random variables 

Consider the entries of {X,} as being independent and identically dis¬ 
tributed random numbers. As we have exchangeable entries, the symmetry 
of the sequence determines th e pr obability of X j being a record to be 1 /j at 
the j th trial; see Renyi (1962) or Arnold et al. (1998). In other words, any 
of the j entries till Xj, including itself, is equally possible to be the record. 
Consider temperatures in a particular city, despite all the fluctuations, it is 
clear that it beats its own record. In formal language, the process of setting 
a record is persistent. Let us denote the number of records in a sequence of 
length n by N n , then the above statement can be rephrased as N n —y oo with 
length of sequence n —* oo. 

The next natural question is about the frequency with which records are 
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broken. Define an indicator function 


J 1 if there is a record at j , 
( 0 otherwise. 


( 1 ) 


Calculating the expected value of Ij is equivalent to calculating the prob¬ 
ability of Ij taking the value 1, but this is precisely the probability of Xj 
being the record- which is 1/j. Thus (Ij) = 1/j. Similarly, for variance, 
we note that expected value of /J is the same as expected value of Ij. This 
immediately gives the variance as 1/j — 1/j 2 . It can be easily proven that 
Ij are pairwise uncorrelated (statistically independent). In other words the 
probability of the position of the records is a Bernoulli process, Ber(l /j). 

As the total number of records N n in a sequence {Xi, X 2 , ■ ■ ■ X n } is 
YTj=\ Iji expectation and variance can be readily calculated: 

n n i 

m = j>i) = £- = "- 

A ’ =1 „. » , ( 2 ) 

v(i v») = E v 'W) = E--Ey 

j= 1 3= 1 J 3 =1 J 


Here H n is the n th harmonic number. A remarkable, well-known fact from the 
theory of records is that for i.i.d. variables these quantities are distribution- 
free, that is independent of the particular underlying distribution p(x), see 
(1998). For example the average number of records (N n ) = 


Arnold et al. 


H n ~ log(?r)+ 7 , where 7 is the Euler constant, is indeed very small compared 
to the length n of the data set; typically records are rare events in independent 
processes. From the above it immediately follows that the variance also grows 
as log(n). Next section on, we will change the notation from N n to A^j. 


1.2. Random walks and records 

Some of the recent works concerns the behavior of records in correlated 



random walks and the few studies on their record statistics is now very briefly 
reviewed. In case of random walks, xt represents the position of a random 
walker at a time step, t. For discrete time steps and jump-lengths drawn 
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from an i.i.d. symmetric distribution, Majumdar and Ziff showed that the 
probability P(M,N) of M records in t — N steps is given by 


P(M, N ) 


^2iV M + 1^ 2 - 2 N+M -1 


M < N + 1 


( 3 ) 


with mean (M) ~ -^=\ /N. Age statistics of record i.e., how long a record 
survives before it is broken, is given by (l) ~ jjjt ^ 


(M) ~ \J (Majumdar and 
Ziff (2008)). Thus in strong contrast to the i.i.d. case, the number of records 


is very large. 

This makes physical sense, as the random walker is indeed trying to go 
somewhere and each time she has genuinely progressed in a particular di¬ 
rection a record is set. Being diffusive, this process is certainly slower than 
ballistic, and in fact grows in the same way as the distance from the origin. 
That said, it is then interesting and somewhat surprising that it retains the 
property of being universal in that it is independent of the exact distribution 
of the jump-lengths. In fact even when distributions have a fat tail (like 


Cauchy distribution), this result has been shown to be valid by Majumdar 
and Ziff ( 2008[ ). 

In the case of continuous time random walks, i.e. when position of a 
random walker is now observed at equal discrete time steps with time interval 
To, the probability P(M, t ) of finding M records within a given time duration 
t, as derived by jSabhapandit (2011), and attains a scaling form given by, 


P(M,t ) ~ (t/T 0 )- Q ' 2 g a (M(t/r 0 )- a/2 ) 

9a{x) = -x~ (1+2/a) L a/2 (x~ 2/a ), 0 < a < 1, 


( 4 ) 


where L a / 2 (x) is a one-sided Levy stable probability distribution function. 
The moments of records, asymptotically go as 


(M u ) ~ (2/ a ) r (^) ( t/ To y a / 2 . 


T(ua/ 2 ) 


( 5 ) 


The mean-age of the record is 


(!) = (t/M) ~ yj) [ln(t/r„) - *(1 - a /2)](i/r 0 )< 1 -“>, (6) 
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where T(x) is di-gamma function, see Sabhapandit (2011). The parameter 
a e (0,1), appearing as exponent of power law decay of the waiting time 
distribution decides the qualitative behavior of record statistics. 


In another generalization due to Wergen et al. (2011]), the record dis¬ 


tribution for a random walk with discrete time steps but asymmetric jump 
distribution, has been calculated for the model governed by x n = x n _i +£ n + c 
with being the jump length with symmetric distribution while c is the con¬ 
stant drift. For Gaussian jump distribution with variance er, in the limit of 
small drift (smallness is compared with respect to c = 0 case) i.e. ^ , 

the mean number of records and record rate (i.e. probability of n th event 
being a record), is given by, 


(. M n (c )) « + — — (narctan(A/n) — y/n), 


n an 


Pn(c) 


— arcta.n(v / n). 


(7) 


7T n a n 


In the large drift limit, P n (c ) approaches a constant value Wergen et al. 


(2011). Later, this problem was exactly solved for arbitrary value of constant 


bias c and a number of interesting results have been obtained such as the 


qualitatively different regions in the space of c, /i by Majumdar et al. (2012) 


where /i is the Levy index for symmetric stable law of jump distribution. 

1.3. Application of record statistics in superconductors and quantum chaos 
The geometric feature of random systems such as size of the largest cluster 
in percolation on a finite lattice of size N, has been shown to follow Gumbel 


distribution in the large -N limit by Bazant (2000). In type-II supercon¬ 


ductors, temperature independence of the magnetic creep rate i.e. rate of 
change of magnetic held in the sample at a constant temperature for a range 
of temperatures has been understood in terms of record dynamics, namely, 
the dynamical properties of the times at which a stochastic fluctuating signal 
(in this case thermal noise due to non-zero temperature) establishes records. 
This puzzling temperature independence of the creep rate, which at non¬ 
zero temperatures has its origin due to thermal fluctuations, has been sorted 
out by showing that the process of vortex penetration into the sample can 
be described in terms of a Poisson process with logarithmic time argument, 


called the log-Poisson process; a result from record dynamics by Oliveira 


et al. (2005). 


6 

























The application of record statistics in case of quantum chaos has been 


treated in Srivastava et al. (2013). It is known that the eigenstate intensi¬ 


ties in fully chaotic systems with no particular symmetries are conjectured 
to behave exactly as these random vectors subject only to a normalization 
constraint. These are also the statistical properties of eigenvectors of the 
classical ensembles of random matrix theory. For chaotic systems, the appli- 


cability of random matrix theory 

Mehta 

(2004) 

Brody et al. 

(1981 

) has been 

well appreciated for long 

Bohigas et al. 

1984 

). As we will discuss in detail in 

Section [3] for the standard map the brea 

dng of the last KAM torus allows for 


diffusion in phase space and diffusion being connected with random walks, 
we expect and indeed observe the average number of records in eigenvectors 
to go as y/N. In fact this is first known instance where K values at which the 
last KAM torus breaks has been seen in a quantity derived from quantum 
mechanical spectrum. 


2. Records in <5-correlated variables 

2.1. Upper records for complex random vectors 

For a correlated sequence, let the probability density for a record variable 
R , at time t be P(R , t). The average record is given by (R) = f d R RP(R, t ). 
Let P(x i ,..., xn ) be the j.p.d.f. of N random variables. The probability that 


the record at time t, is less than R is given by (for example Srivastava et al. 
p0l3b): 

f R 

Q(R,t) = / dx 1 ---dx t P t (xi,...x t ) ( 8 ) 

Jo 

where Pt(xi,... , x t ) — f P(x i, ..., xjy) dx t+ i • ■ • daqv is the marginal j.p.d.f. 
of the first t random variables . It follows that P(R,t ) = d Q(R,t)/ d R. 

Let us specialize to the case of 8 correlated random variables i.e. sum of 
the random entries is a constant. For example, components of normalized 
complex random vectors z n = (n\'tp) } have the j.p.d.f: 


P(zi, z 2 ,..., z N ) = (T(N)/n 



( 9 ) 


This is also the distribution of the components of the eigenvectors of the 
GUE or CUE (Gaussian or Circular unitary ensembles) random matrices, 


for example see the book by Haake (1991). It is the invariant uniform distri¬ 
bution under an arbitrary unitary transformations on the 2N— 1 dimensional 


7 
































sphere. It is the unique (Haar) measure on S 2Ar_1 . The normalization pro¬ 
vides correlation among the components that becomes weak for large N. The 
intensities ay = \zj 2 being the random variables of interest it is more useful 
to define the j.p.d.f. directly in terms of these: 


N 


P(xi ,..., x N -, u) — F(N)8 


Xi. 


u 


( 10 ) 


, i— 1 


where u is an auxiliary quantity, the actual j.p.d.f. corresponding to u = 1. 

The delta constraint also arises in problems such as the “broken-stick” 
one, wherein a stick is broken into a fixed number of pieces but with random 
lengths. The delta correlation is a very weak constraint, and in fact it is 
not hard to conclude on some reflection that the number of records and its 
statistics remains unaffected by this correlation. The way in which a set of 
random numbers can be correlated by the delta function, gives us a hint: 
take i.i.d random normal variables and normalize them, then the resulting 
set of numbers are uniformly distributed on some sphere and hence are delta 
correlated. But the rank order remains the same on a normalization process 
and hence the record statistics remains the same. Note that we have here 
the new variables ay being distributed uniformly on a simplex rather than 
a sphere: however this does not matter for the records as again the act of 
taking a square does not affect the rank order. Nevertheless we show this 
explicitly here, as we derive the probability of the value of the record itself, 
which does differ from the uncorrelated case. Also our analysis generalizes 
in a significant way an earlier analysis of extreme values in such random 


variables by Lakshminarayan et al. (2008). 

Defining 

rR noo 

Q(R,t;u) = / day • • • dx t / P(xi,..., xn;u) dx t+1 ■ ■ ■ dxjy 
Jo Jo 


leads to 


e su Q(R,t]u) du = 


T(N) 


m =0 


-i) r 


m 


~—smR 


(ii) 


( 12 ) 


Using the convolution theorem, and then setting u = 1 in Q(R,t\u) gives 






Hence 


P(R,t ) = ^(-l) m+1 f * 'WjV- 1)(1 -mi2) JV “ 2 0(l - mi?), 

771=1 


(14) 


the probability density that the record is R at time t. Note that P(R, N) is 
the probability density that the maximum value of the entire data set is R, 


which was calculated for the case of random GUE vectors in Lakslnninarayan 
et al. (2008) and therefore P(R,t) here is a generalization. The piecewise 


smooth probability distribution found there has a similar behavior here. 


It was shown in Lakshminarayan et al. (2008) that P(R,N ) is Gumbel 


distributed asymptotically. In fact the generalization presented in Eq. 
is also Gumbel distributed for large N, as for large N and large f > 1 


(13) 


Q(R, £)«(! — exp(— NR)) 1 « exp (—t exp(— NR)). 


(15) 


Since the Gumbel distribution is of the form exp[— exp (—(x — a^/b^} where 
Oat and bpf are the shift and scaling. It follows that for the record statistics the 
relevant parameters are a n = log(t)/IV and b^ = l/N. The shift generalizes 
from \og(N)/N for the maximum, while the scaling remains the same. The 
above form also appears in the limit when the correlations are ignored. 

The average value of the record as a function of time is 


(R(t)) 


1 - 


Q(R , t ) dR 


1 

N 


E (-!r +1 

m= 1 


l 

m 



Ht _ 1 1 

iV “ N 

k =1 


(16) 

where H t is a Harmonic number as defined above. Known asymptotics of the 
Harmonic numbers implies that 


= v (i' +lii w + j • 




( 17 ) 


k =1 


where Bok are Bernoulli numbers, and 7 is the Euler-Mascheroni constant. 
Again, this presents a generalization of the average maximum intensity found 


m 


Lakshminarayan et al. (2008) which corresponds to t = N. 


It is not hard to prove that for intensities of random states too, the 
probability of the position of the records is a Bernoulli process although 
they are correlated by the normalization constraint. Let there be records at 
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positions (ji — 1 < j 2 < ■ ■ ■ < j m ) and let Ij k = 1 if there is a record at jk 
or 0 otherwise. Then the j.p.d.f. 


r m ] 

Prob(/j 1 = 1, I h = 1,..., I jm = 1) = / P(x i, ..., x N ; u = l)dx 1 ■ ■ ■ dx N = TT - 

Jc Z i J 


k =i 
(18) 


1 

jk 


U J3 1 


Here C is the set of constraints: 0 < x^ < Xj 2 , j\ < k < j 2 — 1; 0 < x^ < x 
J 2 < k < j 3 - 1; • • •, 0 < x k < x jm , j m _ i < k < j m - 1; 0 < x k < 1, 
jm < k < N. The above result follows on using the Laplace transform to 


free the constraint in Eq. (10). However this is the result for i.i.d. random 


variables, and implies that the occurrence of a record at jk is an independent 
process, as the above is valid for all arbitrary choices of the locations jk- 
Hence Prob(/j = 1) = 1/j and Prob(/j = 0) = 1 — 1/j, in other words the 
process is Ber(l/j). 

The average number of records is thus 


N 


N 


(N R ) = (j2 i i) = J2i = H «’ 


i j=i 


j= 1 


J 


(19) 


while as a random variable Nr has a distribution essentially given by the 
signless Stirling numbers of the first kind, some times called the Karamata- 
Stirling law (see Nevzorov (1987)). Such laws hold for a variety of disparate 


processes including the number of cycles in a random permutation of N 
objects, number of nodes in extreme side branch of random binary search 
trees etc. Bai et al. (1998). Being distribution-free, the number of records is 


a statistics that directly detects correlations. 

The probability that the final record, which is the maximum in the en¬ 
tire data sequence, lasts for time m can also be simply calculated: denoted 
Sn{jTi) = P{In = 0 ,/ at-i = 0, • • • ,lN-m+ 2 = 0, In-iu+I = 1) = ^/N, it is 
(somewhat surprisingly) independent of m, and uniform. This implies that 
the position at which the maximum occurs is uniformly distributed. The 
implications of this for quantum chaotic wavefunctions where strong scar- 


ring effects of classical periodic orbits ( see 

McDonald and Kaufman 

(1979); 

Heller (1984) 

Biswas and Jain 

(1990); Kudrolli et al. 

(1995) 

Laurent et al. 


(2007)) can affect the maxima of states is of natural interest, and will be 


discussed below. 
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2.2. Lower records in complex random vectors 

The question of a record minimum can be asked in a similar way, the 
cumulative density function for a record minimum, i.e. the probability that 
the record is “greater” than R at a time t is given by 

pOO 

Q(R,t) = / da;i---d x t P t (xi,...,x t ) (20) 

JR 

where P(x i,... , xn ) is j.p.d.f. of N random variables, and Pt(xi,. .., x t ) is 
the j.p.d.f. of t random variables given by j g[[ range P(x i,..., xn) da^+i ... dopy- 
From Q(R,t) we can get P(R,t) as 


P(R,t) 


dQ(R,t) 
d R 


( 21 ) 


ft is not difficult to calculate the probability distribution of record min¬ 
ima in the case of a complex random vector, essentially following the same 
technique as in the case of maximum records, 


Q(R, t; s) = T(N)e~ s ^ ri 
Y{N) e~ sRt 


n 


drje 


i =1 JR 


' N 

n 

i=t +1 ' 


dr t e 


gN-t gt 


Using the convolution theorem and then replacing u = 1, we get 

Q(R , t) = (1 - - Rt) (22) 


Density, 

P(R, t) = — — (N — l)t(l - Rt) N ~ 2 G( 1 - Rt) (23) 

dli 

and average record as a function of t can be easily found 


(R(t)) 


J RP(R,t)dR 
1 

Nt' 


Again this generalizes the results obtained in Lakshminarayan et al. (2008) 


For large N limit Q(R,t ) is again exponential and for t = N, it retrieves all 
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the results obtained in Lakshminarayan et al. (2008). In particular while the 
average intensity is 1/TV, the average minimum intensity is 1/N 2 correspond¬ 
ing to t = N. The record low intensity is inversely proportion to the “time”, 
which in the case of static wavefunctions is the index of the basis used. It is 
easy to show (by just changing the conditions and hence the limits in appro¬ 
priate integrations) that these are also Bernoulli process, Ber (1 jn) and hence 
all the results for survival probability, lifetime distribution etc. remains same 
in this case too. We remark that the exponential form of Q(i?, t) in the large 
N limit is a special case of the Weibull distribution of the standard theory 
of extreme value statistics of i.i.d. random variables. 


3. Record statistics in the classical and quantum standard map 

Attention is now turned to the standard map. This is chosen as the 
standard map is a simply defined dynamical system which has a well-studied 
transition to chaos through the usual route of smooth two-degree of freedom 
Hamiltonian systems. It also has a simple and well-studied quantization 
and allows breaking parity and time-reversal symmetries through quantum 
phases and hence allows for studying GUE, GOE, (or CUE, COE), as well as 
intermediate statistics. Before discussing the record statistics for eigenvectors 
we summarize some aspects of the standard map pertaining to quantization 
and the intensity distribution of its eigenvectors. 

3.1. Elementary aspects of the standard map 

We recall for convenience some well-known basic aspects of the standard 
map. The Hamiltonian of the ^-kicked rotor is given by 

2 00 

H(q,p,t) = V — - ^cos(2vrg) ^ 6(t/T-n). (24) 

n =—oo 

On rescaling Tp —y p and T 2 g —y K, this results in the area-preserving map: 

q n + 1 = (q n +p n ) modi, 

K 

Pn+l = Pn - X- SUl(27T<7 n+ i), (25) 

ZTT 

connecting the position and momentum just after kick n ( q n ,Pn ) to that after 
kickn+1. This is the Chirikov- Taylor map or standard map. The fixed points 
(0, 0), (1/2, 0) of this map are stable and unstable respectively for small K. 
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The periodic property of q endows the phase space with the topology of 
a cylinder. It is also well known that the translational boost in momentum 
by n G Z is a symmetry of the standard map, and this allows one to take 
the modulo 1 for momentum as well and hence endows the topology of the 
two-torus to the phase space. Similarly, it also enjoys a discrete symmetry 
of reflection about the centre of square i.e. p —> (1 — p), q —> (1 — q). In the 
limit K = 0, this map is completely integrable and with the increase in K, 
destruction and creation of invariant surface (i.e. KAM tori ) takes place. 

About K ~ 1, the last rotationally invariant KAM torus breaks, allowing 
global diffusion in the momentum space. If the standard map is unfolded 
to a cylinder it displays normal diffusion in momentum for large enough 
K. When K 5, the classical map is essentially fully chaotic. While for 
very small K the phase space seems to constitute only invariant curves, the 
K = 10 case appears to be completely hyperbolic with no traces of islands. 
At K w 1 a well-known transition takes place when the last rotational KAM 
torus breaks and allows for global momentum transport when viewed on a 
cylindrical phase space. When K < 5 the phase space is clearly mixed with 
large stable islands and chaotic regions coexisting. 

To study the signatures of chaos in quantum mechanical spectrum of such 
systems has been the concern of the subject of “quantum chaos” or “quantum 
chaology”. Naturally quantization of such maps were the first step Casati and 


Ford (1979); Berry et al. (1979). Analogous to Eq. (25), Heisenberg equations 
can be integrated to yield similar equation for the operators. The unitary 
operator connecting states separated by a period of the kick is the quantum 
map. For the standard map on the plane, it is easy to see that corresponding 
unitary operator will be (Chang and Shi (1986); Izrailev ( 1999[) ), 


U = exp 


igT \ f iT 

-^ C0S(2 * 9) CXP 


(26) 


Due to periodicity in momentum and cyclic nature of the position variable, 
the natural phase space will be [0,1) x [0,1) so this map needs to be quan¬ 
tized on torus. The periodicity of both the canonical variables implies that 
the Hilbert space of quantum mechanics is finite dimensional. This dimen¬ 
sionality N is related to a scaled Planck constant as N = 1/h and hence the 
classical limit is the large N limit. A description can be found in the notes 


of Lakslnninarayan (2009); Backer (2003) 
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The quantized standard map in position basis on torus takes the form, 


N—l 


U n n' 


N 


J^exp 


m— 0 


('m + /3 ) 2 .(m + /3) 

-nr ---b 2,7n --— [n — n 


x exp 


N 

' KN 

i -cos 

2vr 


N 

2 ir ( 77 / + a) 
N 


(27) 


(3 and a are the phases which a state acquires along position and momentum 
directions respectively. For periodic boundary conditions, {3 = 0 while for 
anti-periodic, (3 = 1/2. For all other {3 values, time-reversal symmetry is 
broken. A similar role is played by a for parity symmetry. We will focus 
for the case, {3 7 ^ 0 and a 7 ^ 0,1/2 when we can expect that both the 
time-reversal symmetry and parity symmetry are broken and the typical 
eigenstates would be like complex random states. 

To develop more familiarity with the intensities in different eigenfunctions 
of the quantum standard map, we have plotted them in Fig. [TJfor K = 0.3, 0.7 
and 5. It is clear that up to K = 1, intensities are (almost) symmetric about 
p = 1/2 despite a being 0.25. This will have an impact when different aspects 
of record statistics of intensity vectors are considered. 


3.2. Records and anomalous transport in the classical standard map 

The transport of momentum in the standard map has been studied ex¬ 
tensively. It is known that in the strongly chaotic regimes the diffusion might 
be normal or anomalous depending on the absence or presence of very special 
orbits, namely the accelerator modes. Thus momentum diffuses either as nor¬ 
mal Brownian motion or if anomalous then as Levy walks (Chirikov (1979); 


MacKay et al. (|1984 )). Recently a detailed analysis on the role of acceler¬ 


ator modes for anomalous diffusion has been carried out using the so-called 
GALI (Generalized ALignment Index) method (Manos and Robnik (2014])). 
In this section, numerical results for various quantities in the record statistics 
for momentum transport in the standard map is presented. The time series 
being studied are the squared deviation of momentum at n th iteration from 
the initial momentum. As have been already mentioned in the introduction 
that it has been shown in Majumdar and Ziff (2C 


the average record for 
random walk or Levy walk (both), goes as ~ \J At/ tt with time, t while the 
variance of record depends on time linearly i.e. as 2(1 — 2/7T )t. It is then of 
interest to ask if this is the case for the standard map as well. 
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Figure 1: For different K value the Husimi function, \(q,p\'ijj)\ 2 , of represen¬ 
tative states and respective intensity plots are shown where \q,p) is coherent 


state which has been calculated following reference Saraceno (1990). Parity 
seems to be hold good for lower K vkfues despite the value of a being 0.25 
[i.e. maximally broken parity). 


















































For large K , when the standard map is strongly chaotic, the distribution 
of position is uniform in range [0,1). This suggests that in chaotic regime 
one can take the position variable as uncorrelated. Then, the standard map 
reduces to an effective one dimensional stochastic map: 


K 

Pn+i =Pn~ — sin(27rr n ) with r n e Uniform[0,1) 

Z7T 


(28) 


This is precisely the random walk with jump distribution taken from a i.i.d. 
distribution. This clearly brings out not only why the results for random 
walk can be expected to describe the results for record statistics for standard 
map in the hard chaos regime but also that any departure from it is due to 
correlations present in the steps of the deterministic dynamical system. For 
mixed phase regime, when a portion of phase space remains quasi-integrable, 
the time evolution of average records is expected to be sub-diffusive. But, 
any deviation in the large K region can be a signature of accelerator modes. 

Rewriting the Eq. 


25 as 


f n— 1 


(Pn 


Po ) 2 = 2^2 ( ^2sm 2 (27 r q i ) + y] sin(27rgi) sin(27rgj) 


(29) 


, i=0 


The record and the number of records for the quantity (p n — p 0 ) 2 ? namely the 
deviation in momentum square, are studied in this paper. Since, the square 
root of deviation in momentum square, as argued previously, will behave 
like random walk, the exponent in average record as a function of time will 
become 1 for true random walk cases. We have plotted the average record as 
a function of time in Fig. [2] and clearly the long-time behavior is a power law 
oc t a . The averaging is done over 10000 trajectories with initial conditions 
chosen uniformly from [0,1) x [0,1). The left panel corresponds to the case of 
mixed phase space or small K. We see the exponent increasing from values 
smaller than 1 to 1 and for K = 4.05 be as high as 1.33. At this value of 
K there is still a mixed phase space regime, with the stable periodic orbit 
at (0, 0) having just undergone a bifurcation and having become unstable. 
This large value of the exponent may then arise from the presence of period-2 
accelerator modes that have not nearly been studied as much as the usual 
period-1 modes. 

The right panel shows the case for K 5 or hard chaos regime. At K = 


10 there is no evidence of any stable region or accelerator modes (Tomsovic 


and Lakshminarayan (2007)) and it is found that a ~ 1.0. However the same 
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Figure 2: Average record of deviation in momentum square as a function of 
time has been plotted on log-log scale. Clearly the long-time behaviour is 
oc t a . (Left panel) K values and the fitted exponent in the power law, are 
K = 0.5, with a = 0.11 denoted by o, K = 2.00, with a = 0.86 denoted by 
□ , K = 3.2 with a = 1.009 denoted by O and K = 4.05, A with a = 1.33. ft 
is known that K = 4.05, there exist an accelerator mode of period 2. (Right 
panel) The exponent a goes to ~ 1 for K = 10 (denoted by □) i.e. in fully 
chaotic region. A slightly higher value of a for all the K corresponding to 
the regime where accelerator modes are present in the phase space (K = 6.4 
with a = 1.61 denoted by o, K = 12.70,0 with a = 1.495, K = 19.0, A 
with a = 1.3788 and K = 25.25,* with a = 1.2992). The dashed lines are 
the fitting curves. 


figure also shows several other large K values that deviates from this and 
has larger exponents. In all these values there are known accelerator modes. 
As we have averaged over 10000 trajectories with initial momentum drawn 
from uniform distribution in [0,1), a few will satisfy the initial condition 
for accelerator modes which explains why the ( R(t )) for ( p n — p 0 ) 2 is not 
increasing quadratically with time. This shows that simple record statistics 
can be used as a signature to find the anomalous transport region. This 
clearly captures the presence of accelerator modes even of higher periods as 
can be seen in left panel of Fig. [2] 

In Fig. [3j we have plotted the average number of records upto time t. 
Again the averaging is done over the ensemble of 10000 trajectories with 
initial conditions chosen uniformly from [0,1) x [0,1). The left panel again 
corresponds to a mixed phase space regime, while the right to hard chaos. 
Once again the same trends are observed and accelerator modes seem to have 
a strong influence on the number of records set. While in the absence of these 
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Figure 3: Average number of records as function of time has been shown 
for various K values. The averaging is done over 10000 trajectories with 
initial conditions for position drawn uniformly from [0,1) and momentum 
from uniform distribution in [0,1). (Left panel) K values and the fitted 
exponent in the power law, are K = 0.5, with a = 0.22 denoted by o, 
K = 2.00, with a = 0.47 denoted by □, K = 3.2 with a = 0.51 denoted 
by O and K = 4.05, A with a = 0.53. (Right panel) For hard chaos 
regime (K = 10 denoted by □ with exponent 0.52) again the (R n ) varies in 
agreement with the analytical results for random walk. In figure, the various 
K values with the fitted straight lines are, K = 6.4 with exponent 0.60, 
denoted by o, K = 12.70, O with exponent 0.57, K = 19.0, A with exponent 
0.55 and K = 25.25, * with exponent 0.54. The dashed lines are the fitting 
curves.The values quoted here are the best fitted exponents. 


special structures and for the case of hard chaos N R (t ) ~ \/t, which is 


m 


agreement with the analytical results obtained in Majumdar and Ziff (2008) 
for random walks. This implies that the average time for which records lasts 
in the classical standard map’s diffusive regime also increases as a ft as in the 
random walk. It will be interesting to study the maximum and minimum 
intervals over which records last, but we have not pursued that here. In 


the case of random walks this has recently been studied by Godreche et al. 
(2014). 

However it should be pointed out that even if the steps lengths were that 
of a Levy walk the scaling was \ft for random walks, while in the case of 
a deterministic dynamical system, the Levy walk leads to easily discernible 
departures. Thus these Levy walks are evidently correlated ones and the 
correlation plays an important role here. Analytical studies of such walks is 
warranted in this context. 


18 










3.3. Record statistics for standard 7nap eigenvector intensities 

Let us briefly recall the intensity distribution expected from random ma¬ 
trix theory for a complex random state. As the expected value of each in¬ 
tensity component x in an N dimensional space is 1/N, it will be convenient 
to transform to normalized variable y = Nx, and in the limit N —> oo its 
distribution becomes the exponential distribution, 

p(y) = e~ y . (30) 

For N = 2048, 2050, 4096, 4098, this is shown in Fig. [4j Except the tail 
region, it matches very well with the exponential distribution. 



Figure 4: Intensity distribution for eigenvectors of standard map at K = 10, 
a = 0.25, 0 = 0.25 for different N = 2048, 2050, 4096, 4098. 

It is interesting to study the tail region, and the distribution of large 
intensities. Is there a pattern in observing larger intensities (“records”) as 
we go along the index of eigenvectors? What is the probability of occurrence 
of maximum intensity in, say the middle of the eigenfunction? These are some 
questions that we will be addressing comparing with the results obtained for 
random states in Section [2j For large K quantum eigenstates of standard 
map follow the CUE/GUE or COE/GOE results depending on the value of 
the phases a and f3. If ft ^ 0 and a ^ 0,1/2 we can expect that both 
the time-reversal symmetry and parity symmetry is broken and the typical 
eigenstates are expected to be like complex random states. 

The wavefunctions are basis dependent and record statistics will in gen¬ 
eral depend on the space in which the eigenfunctions are represented. For 
small values of K we expect there to be many localized states in the mo¬ 
mentum space while being nearly uniformly distributed in the position, as 
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the dominant invariant classical phase space structures are rotational KAM 
tori. However for large K, position or momentum basis will be equivalent 
statistically. Our analysis below is based on the quantum standard map 
on the torus. It is well known that the eigenstates of the quantum standard 
map on the cylinder are exponentially localized in momentum space |lzrailev 


(1999)) and share phenomenology similar to Anderson localization. Study of 
extremes or records in this case may also be interesting, but is not pursued 
further. 




(a) position representation (b) momentum representation 


Figure 5: The average upper record ( R(t )), from the ensemble of eigenstates 
of the quantum standard map. The parameters used are N = 400 and 
K = 10 (highly chaotic), K = 5 (mostly chaotic), K = 1 (mixed phase 
space), and K = 0.3 (mostly regular). The analytical curve refers to the 
random state result in Eq. (17). In all cases a = /3 = 0.25. 


The average value of the upper record as a function of the index for various 
values of K is plotted in Fig. [5a| This agrees well with the random states 
result in the chaotic region for example when K = 10. However there are 
signihcant deviations in the mixed phase space regime. For example when 
K = 0.3, in the position space most of the records are set up by t/N = 0.5. 
It has been pointed out in Srivastava et al. (2013) that this feature is due 
to the weakly broken parity symmetry. There are signihcant deviations from 
the random state even for K = 5 when the phase space is largely chaotic. 
The momentum space average records in the mixed phase space regime lie 
above the random state result and are not affected so much by the weakly 
broken parity symmetry due to their localization (see Fig. 5b). Thus for 
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small K we see much larger records being set than for the position space. 

A similar picture appears with lower records as well, the results of ran¬ 
dom vector is followed for large K{= 10), while for smaller K the records 
themselves do not vary with t as compared to chaotic cases (see Fig. [6]). As 
is clear that upper record and lower record combined will give the variation 
of the values that function assumes, let’s call it range, therefore in the case of 
standard map for smaller K values (close to integrable regime), say K = 0.3, 
average range of the intensities of eigenvectors is very small as compared to 
large K values (chaotic regime), say K = 10 (See Fig. 5a and Fig. [6j). This 
is consistent with semiclassical analysis that for smaller value of K invariant 
rotation KAM tori are not broken and therefore support wavefunctions with 
smaller range. 



Figure 6: The average lower record (R(t ))lower record) from the ensemble of 
eigenstates of the quantum standard map in position representation. The 
parameters used are N = 400 and K = 10 (highly chaotic), K = 5 (mostly 
chaotic), K = 1 (mixed phase space), and K = 0.3 (mostly regular). The 
analytical curve refers to the random state result in Eq. (24). In all cases 
a = P = 0.25. 


As has been previously discussed, the distribution of the upper (lower) 
record at “time” t is Gumbel (exponential) for large N with appropriate 


shift and scaling. It is shown in Fig. 7a (Fig. 7b) that indeed the upper 
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(lower) record for eigenfunctions of the quantum standard map in the clas¬ 
sically chaotic regime is Gumbel (exponentially)-distributed; also plotted is 
the distribution for the “upper (lower) record” when t = N which refers 

intensity, thus recovering the earlier results of 


to the maximum (minimum) 


Lakshminarayan et al. (2008 




(a) upper record (b) lower record 

Figure 7: The distribution of the upper (lower) records when the index is t for 
eigenfunctions of the quantum standard map with K = 10. After re-scaling 
and a shift, the distributions are of the Gumbel (exponential) type. 


The distribution of the position of the maximum intensity in the position 
representation is shown in Fig. [8aJ where one can see a transition to the uni¬ 
form distribution along with the transition to classical chaos. As has been 
discussed above, in the context of the quantity 5% (m ), the probability that 
the final record, which is also the maximum, lasts for time m is indepen¬ 
dent of m for both i.i.d. and delta correlated variables. Thus the uniform 
distribution of the maximum intensity is consistent with this. 

The sharp peak at the center for small K (here K = 0.3) deserves further 
analysis. In Srivastava et al. (2013), it has been pointed out that for small 
K, when there are many narrow classical resonances, a large fraction of 
eigenfunctions are localized on separatrices and have maximum intensity at 
or very close to q — 1/2, although this is an unstable fixed point (zero 
momentum). For example, with N = 400, K = 0.1 and 0.3, about 75% and 
50% of eigenstates are peaked at q = 1/2. This is because at the prominent 
0/1 resonance corresponding to the fixed points, q — 1/2 is near the turning 
point of orbits on which the eigenstates are localized. 
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As K increases there are more prominent resonances and turning points 
move away from q = 1/2. Eigenstates localize in the resonance interiors, 
thereby the maximum intensity shifts away from q = 1/2. This was the 
qualitative picture put forward in Srivastava et al. (2013) as leading to the 
uniform distribution for maximum intensity for large K. This qualitative 
explanation is also well supported by the distribution of the position of the 
minimum in the position representation. As we expect and indeed observe 
that most of the states have their minimum around stable fixed point, he. 
q — 0 (see Fig. 8 b). The transition in classical behavior from integrable 


(.K = 0 ), to mixed phase space where islands of stability coexist within 
the stochastic sea (intermediate K), to fully-developed chaos is thus well 
captured by the quantity Sn(ui), the distribution of the position of extreme 
intensities. 




(a) upper record (b) lower record 


Figure 8 : The distribution of the position of the last record set, which is also 
the maximum in case of upper record and minimum in case of lower record, 
for eigenfunctions of the standard map with N = 400 and for various values 
of K. In Fig. 8 b inset shows a zoomed in view to clearly see the features of 


S N (m) as function of m for various K values. 


It has been shown above that for random N -dimensional states the num¬ 
ber of intensity records increases as in the i.i.d. case, namely as ~ log N + 7 
both upper and lower records. While we expect to see this for the quan¬ 
tum standard map eigenfunctions in the strongly chaotic, large K regime, 
the mixed and near-integrable regimes show a marked departure that were 
discussed in Srivastava et al. (2013) for the case of upper records. The cor¬ 
relations present for small K lead to results that are similar to those for the 
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random walks, or the classical standard map and results in records increasing 
as a power-law with N. For example it was found that it is almost a pure 
power-law with exponent 0.5 at the critical value K = 0.98, exactly like that 
of random walks. 

4. Summary 

We have shown that time evolution of records for squared deviation of 
momentum of the classical standard map follows a power law. In a heuristic 
manner, it has been argued why in the hard chaos regime this is the same 
as the results for discrete random walks with jump distribution being i.i.d. 
random variables. Hence, any deviation from such a result accounts for 
the correlation present in dynamical system, and is clearly borne out in the 
numerical simulations. The deviation for the large K regime is due to the 
presence of accelerator modes, and result in records that maybe relevant for 
correlated Levy walks. Even, in the mixed phase regime the higher period 
accelerator modes seem to affect the records. 

We have reviewed the results on upper records of intensities of corre¬ 
lated random vectors and derived new ones for lower records. Apart from 
deriving the average record, it has been shown that the probability that a 
record appears at an index j is a Bernoulli process, which is the same as for 
i.i.d. variables. The quantum standard map is a good test bed and has in¬ 
creasingly complex spectrum with the system parameter K. For a quantum 
system with random high-lying states, records’ statistics for random vectors 
applies. Also the study of the position of the last record set in the case of the 
standard map, parametrized by K suggests that beyond a certain value of Ii, 
the eigenvectors do become like random vectors insofar as the records of in¬ 
tensities is concerned. This is also consistent with finding where the number 
of records vs N goes through a transition from a power law to logarithmic, 
as K increases. 

We have essentially begun studying classical and quantum dynamical 
systems using novel measures related to the theory of extremes and records, 
and hope to have shown that this is fruitful and interesting enough to pursue 
further. 
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